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CALCULATION ON IONIZING 
H-ATOM H-ATOM COLLISIONS 
K. Omidvar 
H. Lee Kyle 

ABSTRACT 

The Born approximation is used to calculate the cross sections when an hy- 
drogen atom in an excited state with the principal quantum number n collides 
with an hydrogen atom in the ground state. Calculation has been carried out 
when the first atom is ionized while the second atom may have any final state. 
The cases considered are n = 1, 2, 3, and 4. A closure approximation origi- 
nally used by Lodge is also used here for high impact energy calculations. We 
find disagreement with some of the Lodge's numerical results. 

At sufficiently high impact energies the cross sections are inversely propor- 
tional to the energies. Coefficients of proportionalities for all cases are given. 
When the incident atom is in a high principal quantum number such that the re- 
lation (m/M) K i 2 a Q 2 > n 2 » 1 is satisfied, with K 2 the impact energy in rydberg, 
m and M the electronic and the atomic masses, and n the principal quantum num- 
ber considered, the cross section is given as Q 0 = 4(M/m) [hn 4n 2 - l/2)(77a 0 2 yK 2 . 
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CALCULATION ON IONIZING 
H-ATOM H-ATOM COLLISIONS 

I. INTRODUCTION AND FORMULATION 

In recent years astrophysicists and plasma physicists have expressed in- 
creased interest in the reactions that excited hydrogen atoms undergo when they 
collide with other hydrogen atoms. In this paper 1 we consider the following three 
reactions which lead to the ionization of the excited atom: 


H(rv£) + H( Is) - H + 

+ e + H(ls) , 

(!) 

- H + 

+ e + H(n' 1 ‘) , 

(2) 

- H + 

+ e + H* + e , 

(3) 


where rvt are the principal and the azimuthal quantum numbers of the incident 
atom, and n' V are thosb of the target atom after collision . 

Bates and Grilling 2 have shown that the cross sections for the reactions (1), 
(2) and (3) in the Bora approximation are given by 


Q( iv& - c , 1 s -* n 1 ) 


7TS' 


z>j-r 


|T(tv£m-k, ls-n')| 2 qdq 


(4) 


1 


Here, using atomic units, q is the momentum transfer by the excited atom in 
the center of mass system, m is the absolute value of the magnetic quantum num- 
ber of the excited atom; a m = e/ 2 ^ + 1 ) where €=lifm = 0 and e = 2 if m ^ 0, 
and k is the momentum of the ejected electron. The final state of the target atom 
is indicated by n . For reaction (2) the cross sections reported here include a 
summation over all possible bound excited states of the target atom. When the 
target is also ionized as in reaction (3), n ' = k' , and an additional integration 
with respect to k' is carried out. The initial relative velocity is s, and finally 
the absolute values of the transition matrices are given by 

|T(rv£m“* k, Is - Is ) | = — - 

q 



| T(iv£ m-* k, Is - n' V m') 



e i «‘ r | rv£m><n' V 


m 


• l + r Is > 


( 6 ) 


|T(n^m-k, Is -k')| 



(k |e lQ * r | rvfm^^k' |e l< ** r | ls^ 


( 7 ) 


For evaluation of the matrix elements of exp ( iq • r) in these equations be- 
tween m and k the general formula given by Omidvar and Sullivan** is used, 
while for the squared modulus of the matrix elements of exp (if • r) between Is 


2 



and n' V m', summed over V m', the formula given by Bethe^ is used: 
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( n't 


m e 


iQ*rj 


ls> 
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= 2 8 n 7 q 2 + (qn) 2 J x 

Concerning the limits on q without loss of accuracy we can at high incident 
energies let q max "°°» while to second order in m/M, withm the mass of an electron 
and M the reduced mass of the system, 

AE T m AEl m AE 

^l 1 + ^J ’ << 1 • 0) 

AE being the sum of the excitation energies of the two atoms. For the three reac- 
tions (1), (2) and (3) AE is given in rydberg units by k 2 + 1/n 2 , k 2 + 1/n 2 + 1 - 1/n' 2 , 
and k 2 + l/n 2 +1 + k' 2 , respectively. 

Simplification arises by assuming that q min in (4) is independent of the final 
state of the target atom. It will be shown below that this assumption is justified 
when the incident particle has sufficiently high energy. By putting q min - q n in 
this case, summing the right hand side of (4) with respect to all the excited states 
of the target atom, and making use of the closure property of the hydrogen atom 


~(n ~ 1 ) 2 ^ (qn) 2 ] n ~ 3 
(n + 1) 2 + (qn) 2 ] nt 3 
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wave functions, we find that 


Q 



Is - 




X 


<k |e i «* r | rvCm > 


( 10 ) 


where summation with respect to n ' includes excitation to the discrete, as well 
as to the continuum, states of the target atom. 

Lodge, 5 in order to determine q 0 , has assumed that the largest contribution 
to the cross section arises from the excited states of the target atom for which 
k' = k. This then leads to AE = 2k 2 + 1/n 2 + 1, and through (9) q 0 io determined. 

If, at sufficiently high energies, the contributions from (2) and (3) are added, agree- 
ment is found with the value given by (10) using the above value of q 0 . This can 
be seen from Figures 1-3, in the next section, for the C curve is Just the sum of 
(1) and (10). It should be noted, however, that the integrand in (10) for q small 
behaves as a positive power of q. At sufficiently high incident energies it is 
therefore immaterial what value for q 0 is used as long as q Q « 1. 

In contrast to the case of the charged particle-neutral atom collision, the 
dipole 1/r 7 interaction potential does not exist for neutral atom-neutral atom 
collisions. The long range van der Waals potential behaves as 1/r 7 for distances 
larger than 137 times the atomic radius, andas 1/r 6 for smaller distances, and 
is not important for inelastic collisions. The foregoing discussion implies that 


4 



collisions with large inpact parameters are not important. This will have two 
observable effects. First, few low energy electrons are ejected in neutral- 
neutral collisions in contrast to the charged particle-neutral collisions; second, 
the total cross section at high energies is inversely proportional to the incident 
energy, and the collisional cross section is given by a single parameter. The 
energy distribution of the ejected electrons and these parameters will be given 
in the next section. 

We now derive an analytic expression at high impact energies for the cross 
section when the incident atom is in a highly excited state. Since only close col- 
lisions are important in atom-atom collisions, when n 2 » 1, with n the principal 
quantum number of the incident atom, the interaction of the nucleus of the inci- 
dent atom with the rest of the system can to some extent be neglected, and the 
problem resembles the scattering of an electron by the ground state of the target 
atom. This point has been recognized by a number of workers in the field. 6 

However, the presence of the nucleus of the incident atom has the effect of 
eli m i n at i ng the 1/r 2 potential which exists for inelastic electron-atom collisions. 
We can then treat the problem as the problem of the electron-atom collision 
provided we eliminate contributions from the small momentum transfers which 
corresponds to the long range potential. 

The inelastic cross section in the Bom approximation for scattering of an 
electron by an atom when the momentum transfer is between q, and q and all 

1 x 


5 
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states of the atom are excited is given by 

Q 0 in - ^ £""5 [i-|<o|.‘«-'|o> 

with kj 2 the electron energy in rydberg and <0 1 the initial state of the atom. In 
(11) q x must be lar ger than the minimum momentum transfer for the transitions 
that contribute appreciably to the total cross section. For atomic hydrogen and 
provided Ae/k 2 2 « 1, with Ae the excitation energy of the atom, (11) reduces to 



The second inequality is also the criterion for the validity of the Bora approximation. 

We now fix on an expression for q x . In the impact parameter formulation 
for atom-atom collisions it would be physically reasonable and accurate to ignore 
impact radii greater than the radius of the incident, excited atom. For larger 
radii the interaction between energetic neutral atoms will be negligible. Now q x 
is related to this cut-off radius and an analogous and similarly accurate approxi- 
mation in the momentum change integral is to assume that q x - l/na Q< 

To get the total cross section we must add to (12) the contribution from the 
elastic cross section given by 8 Q 0 e ^ = 7na*/3k* , kj 2 » 1. With these con- 
siderations we find for the total cross section 

. 4rr a 2 

Q. = Qo* + Qo" = -j-j- ten 4"’ - 1/2 ] . 

K 1 
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( 13 ) 


For transitions that contribute appreciably to the total cross section Ac in 
(12) lit*s between 3/4 to 3 rydbergs. We can then combine the condition n 2 » 1, 
the inequalities in (12), and the criterion for the validity of the Born approxima- 
tion in the inequalities 


kj 2 a Q 2 > n 2 >> 1 . (14) 

In terms of the energy of the incident atom (13) and (14) can be combined into 

Q 0 = 4(M/n) ('CnW - 1/2 )(wa # 2 )/fc 2 , 

(m/M) K 2 a Q 2 > n 2 » 1 , ( 15 ) 

with M the mass and K 2 the energy in rydberg of the incident atom, and m the 
mass of the electron. Equation (15) shows that for sufficiently high incident 
energies and high principal quantum number of the projectible, n, the ionization 
cross section increases as n 2 . Examination of Equation (15) for n = 3 and 4 
shows that the results obtained from this equation differ about 2% and 1 % respec- 
tively from the results of the detailed calculation shown in Table 1. 

For low incident energies where the inequalities in (15) are not satisfied but 
still n 2 >> 1, the Bom approximation is not valid and the elastic cross section 
becomes comparable or larger than the inelastic cross section. Derivation of 
an analytic expression in this case does not seem to be easy. It can however be 
said that for sufficiently low energies where the inelastic cross section can be 
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neglected compared to the elastic cross section (cf. Figures 1, 2, 3), the total 
cross section as n increases should approach an upper bound. 

H. RESULTS AND DISCUSSION 

Since the azimuthal angles are cyclic in all the integrations, it is evident 
from (4) that for single ionization a triple, and for double ionization a quadruple, 
numerical integration should be carried out to find the corresponding cross 
sections, the integration with respect to the angle of colatitude of it ' being done 
analytically. These integrations were carried out using the Gaussian quad- 
rature. At impact energies above 100 KeV the sum of reactions (2) and (3) 
differed by at most a few percent from the results given by (10), This later 
equation, with the value for q 0 given by Lodge, was therefore used to obtain 
most of the results at higher energies. This saved much computer time. 

In Figure 1 the sum of the cross sections for reactions (1), (2) and (3), when 
both incident and target atoms are in the ground state, ai e compared with the 
measurements of McClure^ and Wittkower et al. ^ The data of Wittkower et al. 
contains a contribution from charge exchange 

H + H - H + + H" , a6) 

- 1 

Between 3. 15 KeV and 63 KeV McClure determined the cross section for (16) 
and subtracted it out. At higher energies he didn't do this but above 100 KeV the 
contribution from charge exchange is unimportant. ® Whittkower's incident beam 
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consisted only of H(ls) atoms while McClure f s beam contained some excited 
hydrogen atoms. However, McClure determined that their presences increased 
his measurements by at most only a few percent. 

The cross sections for (1), (2) and (3) were previously calculated by Bates 
and Griffing 2 for ground state collisions. Our results for reaction (3) and for 
the sum of the three reactions are about 2% higher than the graphical results of 
these authors. 

Also shown in Figure 1 is the curve designated by C, which is obtained by 
using the results of the closure relationship, Equation (10), for the contribution 
of the excited states of the target atom, and using Lodge f s choice for q 0 . Above 
500 KeV we have estimated reaction (3) by subtracting reaction (2) from the 
closure relation (10). It should be pointed out that neither in this case, nor in 
the cases of n - 2 and 3, do we find agreement with the numerical values of 
Lodge. As explained later on, we attribute this to an error in his calculations. 

In Figures 2 and 3 similar curves with the same definitions are given for 
n = 2 and 3 cases. It is seen from Figures 1 and 2 that the closure relationship 
gives the same results as the explicit integrations for energies larger than 
100 KeV. Below this energy and for n = 2 and 3 it underestimates the total cross 
sections due to an underestimation of reaction (2) in which the target is excited. 
Since the experimental results also agree well with those of the Born approxi- 
mation for energies higher than 100 KeV, it seems that the closure relationship 
provides a useful way of calculating the total cross sections. 
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As discussed in the previous section, due to the lack of the 1/r 2 potential, 
the cross sections fall off at high energies as the inverse of the energy E. In 
Figure 4 the product of the energy and cross section, Ex Q(E), is plotted versus 
E. It clearly is seen that all curves approach their asymptotes 'which are hori- 
zontal lines. Thus at high energies we can write 

Q(n) = C n E 1 ( 17 ) 

In Table I we give C n for the various reactions and n = 1, 2, 3, 4 evaluated at 
1000 KeV. These C n and (17) yield usable cross soctions for impact energies 
greater than 500 KeV. For n greater than 4 use Equation (15). 

Examination of (4) indicates that for the cross section to take the form (17), 
the momentum change integral must become a constant, independent of its lower 
bound q min . We have tested this by putting q min = 0 for all energies of the ejected 
electron and also by use of Lodge’s q 0 for q min and no noticeable variation in the 
cross section occurred at high energies. This insensitivity to q min explains the 
accuracy of the Lodge approximation at high impact energies. 

In order to experimentally investigate the dependence of the total cross 
section on the state n of the projectile, McClure ® has suggested the formula 

Q(n) = Q( Is) n a . (18) 
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Then by (17) 



Using this formula, we find a = 0. 917 for n = 2, a = 0. 834 for n = 3, and 
a = 0. 775 for n = 4. These values are in agreement with the findings of McClure 
which indicate that a must be less than one, although his measurements are for 
the comparatively low energy range of 25-100 KeV. 

For the excited atom in states n = 1, 2, 3 Lodge has previously used reac- 
tions (1) and Equation (10) to calculate the curves c shown in Figures 1-3* His 
results are shown by the dashed lines in Figure 5. The solid lines are our re- 
sults for the same calculation. We attribute the disagreement of the two sets 
of curves to an error in Lodges’ calculations. He indicates that his reaction (1) 
calculation for n = 1 disagrees with that of Bates and Griffing at low energy. 
Lodge attributes this to the fact that he has dropped the second order term in 
q min (cf. Equation (9)). We have checked this and have found that this cannot be 
the cause of the discrepancy. We can only assume that Lodge has some error 
in his calculation the effect of which is more pronounced for n > 1. 

Lodge also repeated the calculation but used a plane wave to represent the 
ejected electrons. Considering that this approximation is not particularly ac- 
curate for impact energies as low as 100 KeV, his results in this case are in 
reasonable agreement with our curves shown in Figure 5. 
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In Figures 6 and 7 are shown the energy distribution of ejected electrons 
when the projectile is a H( 2s) atom with energies of 100 KeV and 500 KeV re- 
spectively. The contributions of the three reactions are shown and the sum is 
compared with the electron distribution in an H( 2s ) -proton collision. In the 
latter case most of the electrons are ejected with low energies and dQ/de for high 
e behaves as lA 2 , with e the energy of the ejected electron. Note the much re- 
duced size of the low energy peak in the problem of atom-atom collisions. For 
the 100 KeV projectile, reaction (I), where the target is not excited, dominates 
the distribution at high e and in this region dQ/de also falls off as 1/e 2 for atom- 
atom collisions. The dominance of reaction (1) is of course much more pronounced 
at lower impact energies. However at higher impact energies reaction (3) be- 
comes most impoitant as Figure 7 indicates. Even in this case reaction (1) is 
the larger in the very high energy tail, but this region contributes little to the 
total cross section. Note that in reaction (3) two free electrons are produced, 
both of these are taken into account in Figures 6 and 7. 

In Tables n and m we present numerical values of the cross sections for 
all three reactions and their sum when the projectile is in the principle quantum 
states n = 1 and n = 2 respectively. Note that in Table HI the 2s and averaged 
2p cross sections are usually about the same size. In fact for all n the cross 
section seems not to vary greatly with the angular momentum i of the excited 
electron. 

All calculations were done on the IBM 360-91 at the Goddard Space Flight 
Center. 
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Table I 


Constants, C if for the asymptotic formula Q. = C^E, 
are given for reactions (1), (2), (3) and the sum. 

Units of C are (KeV 77 a 0 2 ). The C '* are obtained from 
cross sections calculated at 1000 KeV and yield cross section accurate to 
within a few percent at energies greater than 500 KeV. 


Reaction 

(i) 

(2) 

(3) 

Sum (1 + 2 + 3) 

(Is - c, Is - n') 

37.5 

17.6 

71.5 

127 

(2 - c, Is - n') 

50.6 

59.2 

129 

239 

(3 - c, Is ■-* n') 

54.1 

101 

161 

316 

(4 - c, Is - n') 

54.2 

136 

J 

180 

370 












Table n 


Cross sections for the ionization of H(ls) 
when it collides with H(ls) in units of ^a 0 2 ). 


Impact Energy 

Reaction 

KeV 

Log (E) 

a) 

(2) 

(3) 

(1 + 2 + 3) 

1 

0 

0.010 

— 

— 

0.010 

2.5 

0.398 

0.098 

— 

— 

0.098 

5 

0.699 

0. 288 

0.000 

0.000 

0.288 

10 

1.0 

0.499 

0.004 

0.001 

0.504 

15 

1,176 

0.569 

0.017 

0.007 

0.594 

25 

1.398 

0.573 

0.056 

0.043 

0.672 

35 

1.544 

0.531 

0. 091 

0.096 

0. 721 

50.1 

1.7 

0.460 

0.115 

0.185 

0.760 

63.1 

1.8 

0.408 

0.120 

0.241 

0.769 

79.4 

1.9 

0.355 

0.118 

0.285 

0.758 

99.9 

2.0 

0.301 

0.111 

0.309 

0.721 

158.5 

2.2 

0,211 

0. 086 

0.296 

0.594 

199.5 

2.3 

0.172 

0. 074 

0.269 

0.515 

316.2 

2.5 

0.113 

0. 054 

0.200 

0.364 

501 

2.7 

0.074 

0.034 

0.137 

0.244 

1000 

3.0 

0.038 

0.018 

0.071 

0.127 

1,585 

3.2 

0.024 

0.011 

0.046 

0. 081 
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Table m 

Cross sections for the ionization of II (2 -t) in collisions with H(Xs) in units of ("a 
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Figure 1. Cross sections for incident H(ls) on H(1s). Curves I, II, and III are due to the reoctions (1), (2), and (3). Curve C is the 
sum of the cross sections using closure relationship. Equation (10), with qg given by Lodge. The experimental points of McClure 
and Wittkower et al. indicated by squares also contain a contribution from chorge exchange, (16). The open circles indicate meas- 
urements where McClure subtracted out this contribution. Above 100 KeV charge exchange is unimportant. The energy in the ab- 
scissa is that of the incident atom in the laboratory system. 
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Figure 2. The same as in Figure 1, but with the projectile in n = 2 states. 
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in principle quantum states n - 1, 2, 3 and 4. 
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Figure 5. A plot o* total cross section Q versus E 
(10), with qQ given by Lodge. The soli 
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Figure 7. The same as Figure 6 but here 





